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ABSTRACT 

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity 
is reconsidered, with the regular and turbulent magnetic fields treated separately. 
The regular magnetic field strength required to support the gas is consistent with 
independent estimates provided energy equipartition is maintained between turbulence 
and random magnetic fields. Our results indicate that a midplane value of Bq = 4 ^G 
for the regular magnetic field near the Sun leads to more attractive models than 
Bq = 2 fiG. The vertical profiles of both the regular and random magnetic fields 
contain disc and halo components whose parameters we have determined. The layer 
at 1 < \z\ <, 4kpc can be overpressured and an outflow at a speed of about 50kms _1 
may occur there, presumably associated with a Galactic fountain flow, if Bq — 2 fiG. 

We show that hydrostatic equilibrium in a warped disc must produce asymmet- 
ric density distributions in z, in rough agreement with Hi observations in the outer 
Galaxy. This asymmetry may be a useful diagnostic of the details of the warping 
mechanism in the Milky Way and other galaxies. We find indications that gas and 
magnetic field pressures are different above and below the warped midplane in the 
outer Galaxy and quantify the difference in terms of turbulent velocity and/or mag- 
netic field strength. 

Key words: magnetic fields - MHD - ISM: magnetic fields - Galaxy: kinematics 
and dynamics - galaxies: ISM - galaxies: magnetic fields 



1 INTRODUCTION 

Models assuming hydrostatic equilibrium in the Galactic gas 
layer are often used to describe the global state of the inter- 
stellar medium (ISM) (Parker 1966, Badhwar & Stephens 
1977, Bloemen 1987, Boulares & Cox 1990, Lockman & 
Gehman 1991, Kalberla & Kerp 1998). Despite matter out- 
flows in the form of Galactic fountain (Shapiro & Field 1976, 
Bregman 1980, Habe & Ikeuchi 1980) and chimneys (Nor- 
man & Ikeuchi 1989), the assumption of hydrostatic equi- 
librium appears to be viable when considering the average 
state of the ISM over large scales or times. Korpi et al. (1999) 
argue that the warm ISM is on average in a state of equilib- 
rium, although not the hot phase within ~ 1 kpc from the 
midplane. 

Magnetic fields play an important role in the pressure 
balance of the ISM (e.g. Parker 1979, Boulares & Cox 1990). 
Both theory and observations of interstellar magnetic fields 
assert that magnetic energy density should be comparable 
to turbulent and thermal energy densities of the gas (e.g. 
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Zweibel & McKee 1995, Beck et al. 1996), with a direct 
implication for a similar balance in the pressures. 

The Galactic magnetic field has two distinct compo- 
nents, the regular field B on length scales of the order 
1 kpc and the turbulent field b with a typical scale of order 
100 pc. Galactic dynamo theory predicts that the strength 
of the regular field depends on magnetic diffusivity, rota- 
tional shear, the kinetic energy density of the ISM, etc. 
(Ruzmaikin, Sokoloff & Shukurov 1988, Ch. VII; Shukurov 
1998). The turbulent magnetic field strength probably scales 
simply with the kinetic energy density of the turbulent gas 
(Beck et al. 1996, Subramanian 1999). Using the simplest 
form of this scaling in a model of the hydrostatic support 
of the ISM in the solar vicinity allows the pressure sup- 
port from B and b to be separately quantified. In this pa- 
per we discuss implications, for the strength and vertical 
distribution, of both regular and turbulent magnetic fields, 
resulting from hydrostatic equilibrium models. The validity 
of the model is assessed using independent information on 
the magnetic fields in the ISM. 

Descriptions of the state of the interstellar medium tend 
to concentrate on the situation in the solar vicinity. How- 
ever, the gas disc is warped in the outer Galaxy (Binney 



© 2001 RAS 



2 A. Fletcher and A. Shukurov 



1992), and the warping implies that there is an additional 
force involved in the vertical equilibrium balance. In Sec- 
tion ^| we use simple arguments to illustrate some effects of 
the warp on the vertical gas distribution and to suggest that 
the vertical gas distribution can have a characteristic asym- 
metry in the warped disc. Because of the warp, gas densities 
at ±0.5 kpc above and below the displaced midplane can dif- 
fer by a factor of 1.2. An asymmetry of this kind has been 
observed in the outer Galaxy (Diplas & Savage 1991), but 
went unnoticed. 



2 HYDROSTATIC EQUILIBRIUM IN THE 
SOLAR NEIGHBOURHOOD 

2.1 A model of hydrostatic support 

We assume that, on average, the gas at the solar position is 
in a state of vertical hydrostatic equilibrium, 

dP 



Table 1. Gas components of the modelled ISM at the solar po- 
sition 



= pg 



(i) 



where P and p are the total pressure and gas density respec- 
tively, g is the gravitational acceleration and z is the height 
above the midplane. The assumption is probably weakest 
when considering the hot, ionized gas which presumably 
originates in the disc but fills the Galactic halo, and so is in- 
volved in systematic vertical motions. However, hydrostatic 
equilibrium is still a plausible expectation for large areas, of 
order a few kiloparsecs in size. 

The equilibrium model includes several contributions to 
the total pressure, 



Pm ~\~ Pel 



with 



/ J i 



the gas pressure consisting of the turbulent and thermal 
components (with contributions from individual phases of 
the ISM labelled with index i), P m the magnetic pressure 
due to the regular, B, and random, b, magnetic field com- 
ponents, and Per the cosmic ray pressure. Here v Z i is the 
vertical component of the rms turbulent velocity in the i'th 
phase, m the gas number density, T; the temperature and k 
is Boltzmann's constant. 

We separate the pressure due to the regular and tur- 
bulent magnetic fields by assuming that the energy density 
of the turbulent magnetic field is a multiple of the kinetic 
energy density of the turbulent gas, 



B j- b 

8tt 



8 



i -I ^ — > 2 



(2) 



with K a constant and vf — 3v zi is the total turbulent veloc- 
ity. A value of K ~ 1 is in line with theoretical arguments 
showing that equipartition should exist between turbulent 
kinetic energy and magnetic energy densities (Zweibel & 
McKee 1995, Beck et al. 1996). If the turbulent magnetic 
field is concentrated into flux ropes and within these ropes 
energy equipartition also exists (Subramanian 1999), then 
K is the filling factor of the ropes. 

By taking B to be parallel to the Galactic plane we ne- 
glect the contribution of magnetic tension due to the regular 
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Midplane number density no, scale height h, temperature T and 
one-dimensional turbulent velocity dispersion v z for the ISM 
components are given at the solar position. All components are 
assumed to have a Gaussian distribution except the halo gas, 
which is take n to have an expone ntial d istribution. Refe rence s: 
(a) Bloemen ( 1987 ), jbj R eynolds ( |l990a| ), (c) Pietz et al. ( L99j ), 
(d) Kalberla et al. (Ki). 



field to the support of the ISM (cf. Boulares & Cox 1990). 
The vertical regular magnetic field is definitely negligible 
near the midplane away from the disc centre (Ruzmaikin 
et al. 1988, Beck et al. 1996). The dominance of the hori- 
zontal magnetic field is more questionable in the Galactic 
halo. Significant vertical magnetic fields have been detected 
in the haloes of only two galaxies with strong star formation 
(NGC 4631, Hummel et al. 1988, 1991 and M82, Reuter et 
al. 1994), but the Milky Way hardly belongs to this type. 
Vertical dust filaments abundant in many galaxies in the 
disc-halo interface (Sofue 1987) may indicate significant lo- 
calized vertical magnetic fields there. These can be a result 
of shearing the horizontal field in the disc by a galactic foun- 
tain and/or chimneys (cf. Korpi et al. 1999; see also Elstner 
et al. 1995). Still, this interpretation is compatible with a 
predominantly horizontal regular magnetic field in the main 
parts of both the disc and the halo away from the interface 
(e.g. in M51 - Berkhuijsen et al. 1997). In the absence of 
direct observational evidence for the halo of the Milky Way, 
we refer to numerous theoretical models which invariably 
indicate the dominance of the horizontal regular magnetic 
field under a variety of physical setups (e.g. Brandenburg et 
al. 1992, 1993; Brandenburg, Moss & Shukurov 1995). 

The ISM has a multicomponent structure; parameters 
of the basic phases are given in Table Molecular hydro- 
gen, together with the cold and warm neutral phases make 
up the bulk of the gaseous disc mass, the other components 
being the diffuse ionized gas (also known as the Reynolds 
layer), a neutral halo and ionized halo gas (Ferriere 1998, 
Kalberla & Kerp 1998). Different phases do not necessarily 
have a distinct physical nature and the distinction is often 
merely conventional; for example, neutral and ionized halo 
components plausibly represent a single physical entity; like- 
wise, the separation of the warm and diffuse ionized phases 
may simply describe the increase of ionization degree with 
height. The halo components are present in the disc as the 
hot phase. The filling factor for them was taken to be that of 
the hot phase, 0.3 in the midplane, rising to 0.9 at z = 4 kpc 
and 1 further up in the halo, in an approximation to Fig. 10 
of Ferriere (1995) (see also Rosen & Bregman 1995, Korpi 
et al. 1999). To allow for helium, nine percent by number of 
hydrogen atoms has been added to obtain the mass density 
in each phase, and we assume that the helium is fully ionized 
in the ionized halo component. 
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We note that observational estimates of the total ran- 
dom velocity in the diffuse ionized gas yield supersonic val- 
ues 20kms _1 (Reynolds 1990a). This can be due to system- 
atic (albeit random) upward motions of the hot gas forming 
the base of the galactic fountain (cf. Korpi et al. 1999) rather 
than to genuine turbulent motions; the ram pressure of these 
motions contributes to the vertical support of the gas and 
we can include it in the turbulent pressure. 

The pressure of cosmic rays is assumed to be equal to 
that of the total magnetic field, P CI — P m where P m is given 
by Eq. (H). Altogether, the total pressure is given by 



(3K + 1) 



x B 

+ > mkTi + — 



(3) 



2.2 Observational constraints on magnetic field 
near the Sun 

The vertical distribution of B is poorly known, and we derive 
it in Section ^| from Eqs (|l]) and (g), and then assess the 
validity of our results using the following reasonably well 
established constraints. 

First, if equipartition between the energy densities of 
cosmic rays and magnetic fields is assumed then the total 
interstellar magnetic field, B to t = V B 2 + b 2 , at z — is 
estimated from synchrotron intensity as Btot — 6 /iG (Beck 
2001). With observations of synchrotron polarization show- 
ing Bo /Btot — 0.6, as described below, this gives Bo — 4//G 
(Beck 2001); if anisotropic turbulent magnetic fields are 
present this is an upper limit on Bo- Observations of Fara- 
day rotation measures from pulsars and extragalactic radio 
sources yield a midplane value of B of Bo ~ 2 piG (Rand 
& Kulkarni 1989, Lyne & Smith 1989, Rand & Lyne 1994, 
Frick et al. 2001). Note that Bo can be significantly larger at 
a distance of a few kiloparsecs from the Sun because the Sun 
is located close to a reversal of the large-scale magnetic field 
(e.g. Beck et al. 1996). Magnetic field strengths obtained 
from rotation measures can be underestimated if the ISM 
is in pressure balance, so B and p are anti-correlated (Beck 
2001). Therefore we consider models with Bo equal to 4pG 
and 2 /j,G. 

Second, the ratio of the turbulent to regular mag- 
netic field, b/B, is directly related to the degree of polar- 
ization, V, of the Galactic synchrotron emission via V = 
V B 2 /(b 2 + B 2 ) with Vo « 0.7, provided any depolarization 
effects can be neglected (Burn 1966, Sokoloff et al. 1998). 
Berkhuijsen (1971) found b/B ~ 1.2 and Philipps et al. 
(1981) argue that b/B cannot be much larger than unity. Ob- 
servations of the Galactic polarized synchrotron background 
yield 1 <, b/B <, 3 (Spoelstra 1984). This is corroborated by 
analyses of Faraday rotation measures of extragalactic radio 
sources and pulsars which imply b/B ~ 1.7 (Ruzmaikin et 
al. 1988, Ch. IV and references therein) and b/B ~ 3 (Rand 
& Kulkarni 1989) (these are upper limits if B is underesti- 
mated). 

Finally the vertical distribution of the total interstel- 
lar magnetic field, can be obtained from the distribution of 
Galactic synchrotron emission (Phillips et al. 1981, Beuer- 
mann, Kanbach fc B erkhuijsen 1985). The model used by 
Phillips et al. (1981) requires halo emissivity of about 10 
percent of that in the disc, extending to about lOkpc be- 
yond the midplane. Roughly, emissivity is proportional to 
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Figure 1. The regular magnetic field, B, required for precise 
hydrostatic equilibrium with a midplane strength Bq ~ 4 fiG, for 
K = 0.8 (solid), K = 0.6 (dotted) and K = 1.0 (dashed), where 
K is the ratio of turbulent magnetic to kinetic pressures defined 
in Eq. (§). 



the total magnetic field squared or to a higher power if 
detailed equipartition between magnetic fields and cosmic 
rays is maintained (Beck et al. 1996, Zweibel & Heiles 1997, 
Sokoloff et al. 1998). So the total magnetic field, B to t, in the 
halo is expected to be approximately 1/3-1/2 of that in the 
disc, Btot (halo)/B t ot (disc) = 1/3-1/2. 



3 THE ROLE OF MAGNETIC FIELD IN THE 
SUPPORT OF THE GAS LAYER 

The vertical distribution of the interstellar magnetic field is 
even more difficult to obtain from observations than that 
of the gas components of the ISM. Therefore we obtain B 
as a function of z from a hydrostatic equilibrium model 
and then compare the result with independent information 
briefly summarized above. This also allows us to reassess the 
validity of the equilibrium model itself. 

Equation (|l]), with P from Eq. (^), is solved for B(z) 
using different values of K. We use g for the solar vicinity de- 
rived by Ferriere (1998) from Kuijken & Gilmore (1989a,b), 



' 9 = A \z 2 + z zir 2 +M i 



(4) 



where = 4.4x 10" 9 cms" 2 , A 2 = 1.7x 10~ 9 cms" 2 , Z ± = 
0.2 kpc and Z 2 = 1 kpc. 



3.1 Model with B ~ 4/xG 

Constraining the midplane regular magnetic field to Bo — 
4/LtG, as suggested by cosmic ray energy equipartition argu- 
ments (Beck 2001), we obtain the vertical distributions of B 
and b shown in Figs |l| and ^. The best fit to all three con- 
straints described in section Section ^|is with K — 0.8, in line 
with theoretical equipartition arguments (Zweibel & McKee 
1995, Beck et al. 1996); we obtain B = 4/iG, b/B = 1.25 
at z = and B to t (halo)/B to t (disc) ~ 1/3. 

Increasing K amplifies the contribution of the turbulent 
magnetic field and cosmic rays to the equilibrium balance, 
according to Eq. (^). Near the midplane, the strong turbu- 
lent pressure (kinetic and magnetic) created by molecular 
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Figure 2. Turbulent magnetic field strength, 
Fig. (with Bo ~ 4jitG) for if = 0.8 (solid), 
and if = 1.0 (dashed). 



), for the model of 
K = 0.6 (dotted) 



cloud motions is amplified to such an extent that the re- 
quired strength of the regular magnetic field for hydrostatic 
equilibrium begins to drop as K increases (see the solid and 
dashed lines in Fig. |l|) reaching zero for K = 1.4. 

Figure [l] suggests a two-component model of B(z), a 
disc with a scale height of /is ~ 1 kpc and a halo extending 
z ~ 10 kpc with a roughly constant field strength of B ~ 
1 /J.G. The scale height of the diffuse ionized gas (Table |l| is 
roughly the same as Kb- This is also in good agreement with 
Han and Qiao (1994) who have estimated Kb — 1.2 kpc from 
Faraday rotation measures of extragalactic radio sources. 
Polarized radio emission, indicating the presence of a regular 
magnetic field, has been detected high in the haloes of the 
edge-on spiral galaxies to z ~ 4 kpc in NGC 891 and to 
z ~ 8 kpc in NGC 4631 (Hummel, Beck & Dahlem 1991). 

Figure |^ shows a similar two-component structure of 
the turbulent magnetic field 6; a disc with scale height of 
hb — 0.5 kpc and a slowly varying field in an extended 
halo (with b decreasing by a factor of two at z ~ 7 kpc), 
where B ~ b. We note a significant difference between the 
scale heights of the regular and turbulent magnetic fields, 
hs/hb — 2 for the narrower components. This may have im- 
portant implications for theories of galactic magnetic fields. 



3.2 Model with B ~ 2/j.G 

If Bo is significantly smaller than 4 fiG then the best fit 
to the constraints described in Section |2| is with K — 1.6 
(Figs |^ and ^); we obtain B — 2.4 /^G at the midplane 
(but exclu ding the molecular gas; see the discussion in Sec- 
tion p^ ), b/B w 2.4 and fltotal (halo) /B tota i (disc) « 1/3. 
Realistic variations of the gas parameters do not produce 
outcomes meeting all three criteria far from K = 1.6, which 
results in a constraint 1.5 ^ K ^ 1.7. Using a different rep- 
resentation of g for the stellar disc derived by Rohlfs (1977, 
p. 35), with gravity from a spherical dark matter halo added 
as Kot 2 / Rq, where V ro t is the rotation speed of the disc and 
Rq the solar radius, the range becomes 1.3 K 2.4. If K 
is beyond this range, then either Bo > 3 piG, or Bo < 1 /xG, 
or bo I Bo > 3 must be accepted for the midplane values. 
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Figure 3. The regular magnetic field, B, required for precise hy- 
drostatic equilibrium with a midplane strength Bo — 2 /j,G, for 
different values of K with both the pressure and mass contri- 
butions from molecular clouds omitted (solid line) and with the 
molecular clouds included (dotted line up to z ~ 0.2 kpc, solid 
line above). In those ranges of z where B = 0, the hydrostatic 
equilibrium equations formally require negative magnetic pres- 
sure thereby indicating an overpressured region. 



3. 2. 1 Overpressure above z ~ 1 kpc 

A strange feature of Fig. for the better fits at K = 1.4 
and 1.6, is the apparent absence, or a strong decrease, of 
B at z = 1-6 kpc. The reason for this is clear from Fig. ^[ 
where we show the total pressure excess over the hydrostatic 
equilibrium value, 



AP = P 



pgdz 



(5) 



assuming B(z) is as in Fig. |^ but not less than 1 fiG at any 
z. We see that there is more than enough pressure from the 
gas, turbulent magnetic field and associated cosmic rays to 
support the weight of the overlying ISM at 1 ^ z ^ 6 kpc. 
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Figure 4. Turbulent magnetic field strength, b for the model with 
Bo ~ 2 fiG, excluding the effect of molecular cloud turbulent 
energy, for K = 1.2 (dashed), K = 1.4 (solid) and K = 1.6 
(dotted). 
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Figure 5. Pressure excess over the hydrostatic equilibrium value, 
described by Eqs. (^|) and (g), using B shown in Fig. ^ but with 
minimum of B enforced to be 1 fiG, for K = 1.4 (solid) and 
if = 1.6 (dashed). Positive AP indicates overpressure. 



Any further contribution from the regular magnetic field 
enhances overpressure in this range of z. 

The excess pressure at z = 1-4 kpc, if real, appar- 
ently contributes to driving the Galactic fountain flow. 
The vertical velocity driven by the overpressure is V z — 
(|AP/p) 1//2 ~ 50kms _1 , significantly smaller than the 
magnetosonic speed (100 km s -1 or more) and compatible 
with what can be expected for galactic fountains. It is no- 
table that a similar systematic outflow of the hot gas occurs 
in models of multiphase ISM driven by supernovae (Rosen 
& Bregman 1995, Korpi et al. 1999, Shukurov 1999). 



3.2.2 Molecular clouds 

An unsatisfactory feature of the distribution shown with the 
dotted line in Fig. ^ (i.e. for a model including the molecular 
gas) is the deep minimum of B at z — 0. This feature also 
occurs in the model of Boulares & Cox (1990, their Fig. 2). It 
is very unlikely that B can have such a sharp minimum even 
if it has a global odd parity in z, resulting in the constraint 
Bo = B(0) = 0. Apart from increasing Bo to 4/iG (Sec- 
tion 3.1), a possible resolution of this difficulty assumes that 
the coupling between kinetic energy density and turbulent 
magnetic fields, parameterized by K, is weaker for the bulk 



motion of molecular clouds than for other ISM phases. This 
assumption is physically compelling because the clouds rep- 
resent the densest, self-gravitating part of the ISM, rather 
isolated from the other diffuse phases. The minimum in B 
does not occur if K 5s 0.5 for the molecular clouds and 
K ~ 1 in the other phases. If, otherwise, K ~ 1 for the bulk 
motion of molecular clouds as well, then either B must grow 
sharply in the range \z\ < 100 pc from a small midplane 
value, or B > 3 fiG at z — (cf. the upper panel of Fig. ^) . 

A possible interpretation is that the turbulent magnetic 
field interacts differently with the molecular clouds than 
with the diffuse phases of the ISM (e.g. because of the fluctu- 
ation dynamo action in the diffuse gas - see Sect. 4.1 in Beck 
et al. 1996 and Subramanian 1999). If the fractional volume 
swept by the molecular clouds over the correlation time of b 
is significantly less than unity, only a correspondingly small 
fraction of the magnetic field will be coupled to the clouds. 
The number of clouds in a correlation volume I 3 of the turbu- 
lent magnetic field is TV ~ (l/l c ) 3 fc, where f c is the volume 
filling factor of the clouds and Z c is the cloud size. During the 
correlation time of the magnetic field, r ~ l/v, the clouds 
moving at a speed v c w 8kms _1 (cf. Table [j]) sweep the 
volume V ~ NlcV c T. So the clouds affect magnetic fields 
within the fractional volume V/l 3 ~ f c v c l/{vl c ), where v is 
the r.m.s. turbulent velocity in the ambient medium. With 
f c = 10 -3 -10 -4 (Berkhuijsen 1999 and references therein), 
l/l c ~ 10 and v/v c ~ 2, we obtain V/l 3 ~ 5 x (10~ 2 - 
10 -3 ) <C 1, so the coupling between molecular clouds and 
the intercloud magnetic fields is expected to be weak. Molec- 
ular clouds move with respect to the surrounding gas, so 
their connection with the external magnetic field can be 
lost rapidly. Equipartition can still be maintained within 
the clouds but may not occur between the bulk motion of 
the clouds (that provides the average pressure support to 
the ISM) and the turbulent magnetic field. 

The same arguments apply to the regular magnetic field 
and one might expect a relatively weak dependence of B on 
the parameters of the molecular gas (cf. Beck 1991). We 
show in Fig. ^ the vertical profile of the turbulent magnetic 
field b calculated from Eq. (0) for a choice of K values, where 
the contribution of molecular clouds to kinetic energy den- 
sity has been suppressed. With the contribution of molecular 
gas to b suppressed in Eq. (||), the turbulent magnetic field is 
reduced near z = giving room for a stronger regular mag- 
netic field. The scale height hi, is not sensitive to different 
values of K. 



3.3 Equilibrium near the Sun: discussion 

The comparative simplicity of the results obtained with 



Bo — 4/iG, compared to those discussed in Section 3.2 
for weaker Bo, as well as their quantitative and qualita- 
tive plausibility leads us to favour the arguments suggesting 
Bo — 4 /iG in the solar vicinity. 

The results of our model with Bo — 4 ^iG, discussed 
in Section JO], compare reasonably well with the vertical 
distribution of the magnetic field derived by Kalberla & 
Kerp (1998, K&K hereafter) using a hydrostatic equilibrium 
model, shown in their Fig. 5. In particular, K&K predict a 
regular magnetic field in the halo of B ~ 1 pG at z ~ 8 kpc 
as a requirement for hydrostatic equilibrium. However, K&K 
have a weaker turbulent magnetic field, with feo — 3.5 pG, 
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Figure 6. Expected percentage polarization of non-thermal radio 
emission, using the regular and turbulent magnetic field distribu- 
tions from Figs 0and| (for K = 0.8) in Eq. M, with Vo = 70%. 

that extends to only z ~ 2 kpc compared to z > 8 kpc in 
our model. Polarized radio emission from the haloes of the 
edge-on spiral galaxies NGC 891 and NGC 4631 is signifi- 
cantly depolarized by turbulent magnetic fields (Hummel et 
al. 1991) and any large scale fields present in the halo must 
become tangled by turbulence and Parker and thermal in- 
stabilities (Beck et al. 1996), so we expect to see a turbulent 
component of the magnetic field at all z. 

In contrast to K&K, who have magnetic pressure in the 
disc of 1/3 the total gas pressure and no pressure contri- 
bution from cosmic rays, we find equality between the gas 
(turbulent plus thermal) and total magnetic field pressures 
in the region \z\ < lkpc. And whereas K&K have pressure 
equilibrium between gas and magnetic fields in the halo at 
2 ~ 5 kpc our model only requires P m ~ 0.5P g . We have as- 
sumed equipartition between cosmic ray and magnetic pres- 
sures at all 2. 

3. 3. 1 Implications for external galaxies 

The observed degree of linear polarization of synchrotron 
radio emission is given by Burn (1966) as 

p= *ww' (6) 

where B± and b± are, respectively, the regular and turbulent 
(isotropic) magnetic field components transverse to the line 
of sight and Vo ~ 70% is the degree of polarization in a 
purely regular field. Assuming that b\jB\ ~ b 2 /B 2 Fig. ^ 
shows the expected variation in V with z. 

An implicit assumption in Fig. |^ is that Vo is constant 
in 2. However, Vo — (7 + l)/("y+ |) depends on the electron 
energy spectral index 7. If 7 increases with 2, due to cosmic 
ray electron energy losses as they travel into the halo, then 
Vo and hence V will increase with 2 making the slope at 
2 > 3 kpc in Fig. | even steeper, although this effect is only 
weak. 

Observations of polarized radio emission from edge-on 
galaxies, such as NGC 891 (Hummel et al. 1991, Dumke et 
al. 1995) and NGC 4631 (Hummel et al. 1991, Golla & Hum- 
mel 1994, Dumke et al. 1995) show the degree of polariza- 
tion rising with 2 from a minimum at 2 = 0. The observed 
small V at the midplane, e.g. V\ z =o ^ 1% in NGC 891 
(Hummel et al. 1991), can be explained by strong Faraday 



depolarization by both regular and turbulent fields in the 
disc. Above 2 ~ 2 kpc, the steady increase in V shown in 
Fig. ^| is reproduced in the observations, but V in Fig. ^ 
is a factor of two higher than observed in external galax- 
ies at the wavelength A = 20 cm. Faraday depolarization 
due to both regular and turbulent magnetic fields is ex- 
pected to be weak in galactic halos where Faraday rotation 
measure does not exceed lOradm -2 due to low gas den- 
sity. For example, depolarization by internal Faraday disper- 
sion arising from a turbulent magnetic field b ~ 2 fj,G with 
correlation length I ~ 0.5 kpc (Poezd, Shukurov & Sokoloff 
1993), path length through the halo L ~ 10 kpc and ther- 
mal electron density n c — 10 -3 cm -3 can only reduce the 
degree of polarization by a few percent even at A = 20 cm, 
V = Vo(l - e- s )/S w 0.96P , where S = 4A 4 (0.81n c fo) 2 /L 
(Burn 1966, Sokoloff et al. 1998). 

A plausible explanation of the relatively low polariza- 
tion of synchrotron emission from galactic halos of edge-on 
galaxies is that b\jB\ is larger (~ 6) than predicted by our 
model (2-3) for the Milky Way. 



4 HYDROSTATIC EQUILIBRIUM IN THE 
OUTER GALAXY 

The Galactic gas disc is warped beyond the solar orbit (Bin- 
ney & Tremaine 1987, Binney 1992). The warp is thought 
to propagate through the disc in azimuth and the speed 
of propagation determines whether the ISM at large radii 
has enough time to reach a state of hydrostatic equilibrium. 
The period of the warping wave is about 3 x 10 8 yr, assumed 
to be similar to the disc rotation period at a radius 12 kpc 
where the warp starts (Binney 1992). This is longer than the 
sound crossing time over a scale height of 1 kpc at the sound 
speed of lOkrns -1 . Thus, hydrostatic equilibrium can be a 
valid first approximation for the state of the gas at R = 12- 
18 kpc (see, however, Masset & Tagger 1996, for a discussion 
of possible deviations from equilibrium). 

In this section we argue that the warping can lead to 
a noticeable asymmetry in the vertical profile of the gas 
density and confirm this using the vertical H I distributions 
in the outer Galaxy obtained by Diplas & Savage (1991). 

We describe positions in Galactocentric polar coordi- 
nates (R, 0,2), where R is distance in the Galactic plane 
from the Galactic centre, O is the azimuth in the Galactic 
plane, with the Sun at O = 180° and increasing anti- 
clockwise viewed from the north Galactic pole; and 2 is ver- 
tical distance from the Galactic plane. The galactocentric 
distance of the Sun is adopted as _Rq = 10 kpc, rather than 
the currently used Rq — 8.5 kpc to avoid the need to rescale 
the data of Diplas & Savage (1991). 

4.1 The observed vertical asymmetry in Hi 
density 

Data from the northern hemisphere Hi survey (Stark et al. 
1991) was used by Diplas & Savage (1991) to study the gas 
morphology in the outer Galaxy. The data show a distinct 
asymmetry in the vertical distribution of Hi across a sub- 
stantial region of the outer disc. The asymmetry can also 
be seen, but less distinctly, in maps of the outer Galaxy 
presented by Burton (1988, Fig. 7.18). 
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Figure 7. Hi number density from Diplas & Savage (1991) at 
various R and © = 90° (diamonds) and fits using two different 
velocity dispersions in a warped gas disc with the asymmetric 
gravity of Fig. ^. Dashed: u = 7, 8 and 5kms _1 from top to 
bottom. Dotted: and u = 15, 15 and 8 kms" 1 from top to bottom. 



The asymmetry occurs over a sector of the outer Galac- 
tic disc in the azimuthal range B = 60° to G = 110°, visible 
between R — 14 and 18 kpc. It is notable that the warping 
is maximum in this region (Diplas & Savage 1991). Unfor- 
tunately, the region around O = 270° where the other crest 
of the warp is located is not accessible for a northern hemi- 
sphere survey. 

Diplas & Savage (1991) do not comment directly on the 
asymmetry, but the data reduction process used can lead 
to some mis-positioning of gas, especially at lower densities, 
because of confusion with Hi clouds moving with peculiar 
velocities. (B. Savage, 1998, private communication). Re- 
versing the calculations of Diplas & Savage (1991) in the 
region of asymmetry shows that the velocity of the Hi , 
— 100 < v < — 40kms _1 , is indeed similar to that of in- 
termediate velocity clouds (Wesselius & Fejes 1973, Kuntz 



Figure 8. The vertical acceleration in a warped disc. Radial 
extrapolation, to R = 16 kpc, of the acceleration due to grav- 
ity from the solar position after Ferriere (1998), \g\ from Eq. (^) 
(dotted) and \g + Qol with constant ag added to move the mini- 
mum to z = 0.5 kpc in a rough model of the force responsible for 
the warp (solid). 



& Danly 1996). However the region (35° ^ / ^ 78° and 
—8° b ^ 8°, where I and b are Galactic longitude and 
latitude) is too large to be associated with a single cloud. 
We proceed with the assumption that the asymmetry is a 
real feature of the outer Galaxy gas distribution. 

The results of Diplas & Savage (1991), for three po- 
sitions in the outer Galaxy, are illustrated in Fig. The 
maximum gas density is displaced to positive z, reflecting 
the warp of the gas disc. It is also clear that the gas density 
is asymmetric with respect to the displaced midplane, with 
a density excess at positive z. 



4.2 Asymmetry in the vertical H I density caused 
by the warp 

To model gas equilibrium in a warped disc, we describe the 
warp locally as being produced by a ^-independent acceler- 
ation ao, added across all z as required to make z — 0.5 kpc 
the bottom of the potential well. Since the stellar disc is 
not significantly warped, we scaled g of Eq. (JiJ) to the outer 
Galaxy (Ferriere 1998) (neglecting edge effects and gravity 
from the gas since the region of interest is not far from the 
edge of the stellar disc) and then added a constant accelera- 
tion ao = 1.5 x 10 -9 cms -2 at R = 16 kpc. The modulus of 
the resulting vertical acceleration is shown in Fig. |§| Acceler- 
ation of this magnitude displaces the gas density maximum 
to the observed height at this radius. Thus, the vertical equi- 
librium equation (hj) is replaced by 

dP 

— = p(g + a ) , 
with 



(7) 



-9 : 



Aiz 



exp 



R - R r, 

Ri 



z 



R% + Ro 



R 2 



R\ 



(8) 



and 



where A±, A2, Z\ and Z2 are defined after Eq. 
Ri = 4.9 kpc, R 2 = 2.2 kpc. 

If all of the pressure sources are symmetric with respect 
to the shifted midplane, the asymmetric force of Eq. (Q) 
should give rise to a larger gas scale height above the shifted 
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midplane at 0.5 kpc, similar to what is observed - see fits in 
Fig. 

To further, investigate the effect of the warp, we use 
Parker's (f966) model of hydrostatic equilibrium 



p(z) = po exp 



-*(*) 



(1 + a + /3)u 2 



(9) 



where 4> is the gravitational potential, obtained by integrat- 
ing Eq. (^) and a and fi represent the pressures due to mag- 
netic fields and cosmic rays, respectively, as multiples of the 
gas pressure. A more detailed analysis as in Section ^ is 
hardly warranted here because of our limited knowledge of 
the gas parameters in the outer Galaxy. Assuming a = /3 = 1 
and using the asymmetric force of Eq. (j?]) , we fitted Eq. ^ 
to the observed density distribution using different values of 
the effective total gas velocity dispersion, u. 

As shown in Fig. [?], the density distributions are re- 
produced satisfactorily close to the shifted midplane with a 
single (smaller) velocity dispersion. However, there is addi- 
tional asymmetry at larger heights which can be accounted 
for by larger values of u at negative z, at least at R — 14 
and 16 kpc. The same effect is found at other positions in 
the sector of the Galactic disc that is positively warped; at 
each position a further source of asymmetry is required as a 
higher order effect at large heights. 

To quantify deviations from Eq. (Q) with a constant ao, 
we show in Fig. ^| the residual acceleration 

ai = p — g - a , 

oz 

at R = 16 kpc both for a uniform velocity dispersion 
u — 11.5 kms -1 , a mean between the two approximations of 
Fig. ^, and for u = 8 kms -1 (shown dotted in Fig. 0) above 
the shifted midplane and u = 15 kms -1 (dashed in Fig. ^) 
below the midplane. The residual, asymmetric component of 
the acceleration is |ai| ~ 10 -9 cms -2 or less at Az = ±lkpc 
from the midplane of the warped disc, so lai/ao| < 0.7. 

The former option, illustrated in Fig. Ha, attributes the 
whole residual asymmetry to a ^-dependent component of 
the warping force, ai, which then must vary at a scale of 
order 1 kpc, with oi possibly becoming comparable to ao at 
large heights. 

In the latter case, illustrated in Fig. ^b, one supposes 
that a part of the residual asymmetry is due to an asym- 
metric pressure in the ISM. We have quantified the residual 
in terms of velocity dispersions. An obvious candidate for 
the residual asymmetry can be the regular magnetic field 
B(z). We note that some galactic dynamo models are com- 
patible with strongly asymmetric B(z) where B(z) vanishes 
in a layer at a height of about 1 kpc on one side of the disc. 
This could happen if regular magnetic fields in the disc and 
halo have different parities in z (Sokoloff & Shukurov 1990, 
Brandenburg et al. 1992). The range of ai in Fig ^|b, Aai, 
would correspond to an asymmetric component of B(z) of 
order \/~pKa\Kz ~ 1 /iG, where Aoi — ai, gas density cor- 
responds to 0.03cm -3 , and B has to be smaller above the 
midplane z > 0. 



5 SUMMARY 

As discussed in Section]^, parameters of a multiphase ISM in 
the solar vicinity are compatible with a model of hydrostatic 
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Figure 9. Solid: \p~ 1 dP/dz\, the total vertical acceleration due 
to gravity and the warping force required to maintain the observed 
H I in hydrostatic equilibrium at R = 16 kpc. Dotted: the accel- 
eration \g + ao\ as in Fig. ^| Dashed: oi, the difference between 
the above two. (a): the above with a single velocity dispersion, 
u = 11.5 kms , and (b): a composite best fit with u = 8kms — 1 
above the shifted midplane and u = 15 kms -1 below. 



equilibrium where approximate equipartition exists between 
the energy density of the turbulent magnetic field and the 
turbulent energy density of the gas in its diffuse phases. 
Models of hydrostatic equilibrium are sensitive to the mid- 
plane strength of the regular magnetic field, Bo- 

If Bo ~ 2 pG, the ISM may be overpressured at 1 <, z <, 
4 kpc (so hydrostatic equilibrium is a poor approximation 
there), and this can contribute to driving the galactic foun- 
tain. The interstellar magnetic field should be coupled more 
weakly to the molecular clouds than with the other, diffuse 
gas phases; otherwise, the regular magnetic field strength 
has a deep, narrow minimum at z — 0. 

If Bo ~ 4pG then our model gives two-component ver- 
tical profiles of both regular, B, and turbulent, 6, mag- 
netic fields. The narrower components have scale heights 
of Jib — lkpc for the regular magnetic field, hi, ~ 0.5 kpc 
for the turbulent magnetic field. The more extended compo- 
nents have an approximately constant regular magnetic field 
strength of B ~ 1 fiG for 2 <, z <, 8 kpc and bo /Bo = 1.25 
at 2 = 0. The vertical dependence of the degree of linear 
polarization of synchrotron emission expected from the de- 
rived B(z) and b(z), shows a steady increase above z ~ 2 kpc 
similar to that observed for external edge-on galaxies. 

Models with Bo — 4 /xG seem to be simpler and more 
attractive than those with Bo — 2 /jG. The former are consis- 
tent with hydrostatic equilibrium maintained at all heights, 
whereas the overpressure at 1 ^ z ^ 4 kpc in the latter mod- 
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els drives systematic vertical gas flow to the halo at a speed 
of about 50kms _1 . This counterintuitive conclusion results 
from a larger contribution of turbulent magnetic fields and 
cosmic rays to the total pressure required in the latter mod- 
els for hydrostatic equilibrium at z Si 1 kpc. We should em- 
phasize that our results with So — 2 fiG are only tentative 
at z <; 1 kpc where our basic assumption of hydrostatic equi- 
librium does not apply with this value of Bq. 

We have demonstrated in Section ^ that the vertical 
force acting on the gas in a warped disc should become asym- 
metric in z, resulting in asymmetry in the gas distribution. 
By taking the simplest form for the vertical structure of the 
warping force as a first approximation and assuming hy- 
drostatic equilibrium can be reached, we calculated the ex- 
pected vertical gas density distribution. The observed warp- 
ing is produced by an acceleration of do = 1.5 x 10 -9 cms -2 
at R — 16 kpc. The vertical profiles of neutral hydrogen den- 
sity in the outer Galaxy (Diplas & Savage 1991) show an 
asymmetry of the expected form, with the density 0.5 kpc 
above the shifted midplane a factor of 1.2 higher than the 
density 0.5 kpc below. 
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